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Abstract 

We study the integrable structure of the 2D Laplacian growth problem with 
zero surface tension in an infinite channel with periodic boundary conditions in 
the transverse direction. Similar to the Laplacian growth in radial geometry, this 
problem can be embedded into the 2D Toda lattice hierarchy in the zero dispersion 
limit. However, the relevant solution to the hierarchy is different. We characterize 
this solution by the string equations and construct the corresponding dispersionless 
tau-function. This tau-function is shown to coincide with the genus-zero part of 
the generating function for double Hurwitz numbers. 
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1 Introduction 



Growth problems of Laplacian type such as Hele-Shaw viscous flows refer to dynamics of a 
moving front (an interface) between two distinct phases driven by a harmonic scalar field. 
These essentially nonlinear and non-local problems attract much attention for quite a long 
time pp. The Laplacian growth problem appears in different physical and mathematical 
contexts and has important practical applications. The most known ones are filtration 
processes in porous media, viscous fingering in the Hele-Shaw cell, electrodeposition and 
solidification in undercooled liquids. The problem in 2D is the most studied one. For 
reviews see [21 E]- To be definite, we shall speak about the dynamics of an interface 
between two incompressible fluids with very different viscosities. In practice the 2D 
geometry is realized in the Hele-Shaw cell - a narrow gap between two parallel plates 
(Fig. I}. 

Remarkably, 2D Laplacian growth (LG) with vanishing surface tension possesses a 
hidden integrable structure which, for the problem in the radial geometry, i.e., in the 
plane with a point-like source or sink, was revealed in (3] and further studied in [5]- jllj . 
Since evolution of planar simply-connected domains is most naturally described by time- 
dependent conformal maps, there is no surprise that this structure is actually immanent 
for general conformal maps and classical boundary value problems. Specifically, it was 
shown in |H [6] that evolution of conformal maps is governed by an integrable hierarchy 
of nonlinear partial differential equations which is a zero dispersion version [121 03] of the 
2D Toda hierarchy [13] . We call it the dispersionless Toda (dToda) hierarchy. The times 
of the hierarchy are harmonic moments of the evolving domain and the Lax function is 
identified with the conformal map. In fact the dispersionless Lax equations for it can 
be derived in the framework of the classical theory of conformal maps depending on 
parameters. 

It is also remarkable that the integrable structure unites the LG problem with im- 
portant areas of mathematics and theoretical physics such as inverse potential problem, 
quadrature domains, random matrices, theory of solitons and c = 1 string theory. Some 
of these links are reviewed in [T5]. We are going to add a new item to the list. 

One aim of this work is to make explicit the integrable structure of the LG problem in 
a different geometry - namely in an infinite channel with periodic boundary conditions 
in the transverse direction (an infinite cylinder). In this version, the problem is also 
known as the Saffman- Taylor problem. We show that the evolution is governed by the 
same dispersionless Lax equations but the specific solution of the dToda hierarchy is 
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substantially different from the one relevant to the LG growth process in the plane. 

Another aim is to show that the solution to the LG in a cylinder has an intriguing 
combinatorial and algebro-geometrical meaning which is more transparent than in other 
settings. Namely, it appears to be closely related to the enumerative geometry of ramified 
coverings of the Riemann sphere. The corresponding tau-function turns out to be the 
dispersionless limit of the tau-function for Hurwitz numbers which is the generating 
function for numbers counting ramified coverings of CP 1 of a certain ramification type (see 
[19] for a review of the Hurwitz problem and related topics). The combinatorial theory 
of ramified coverings of CP 1 was linked to integrable systems in [201 12T] . In particular, 
the generating function for double Hurwitz numbers was shown to be a special solution 
(tau-function) of the 2D Toda lattice hierarchy in [2T]. The integrability of Hurwitz 
partition functions and their relation to matrix models is now actively investigated (see 

e.g. m-m)- 

In the rest part of the introduction we outline the contents of the paper. 

The LG problem on the surface of a cylinder. Consider an infinite cylinder of 
radius R obtained from the physical (X, Y) plane by identifying the points (X, Y+2nmR) 
for all m G Z. As usual, we will use the complex coordinates Z = X + iY, Z = X — iY. 
Let T be a closed non-intersecting contour on the cylinder equivalent to the non-trivial 
cycle. It divides the cylinder into two infinite pieces, D_ and D + , - to the left and to 
the right of V respectively. Without loss of generality we assume that the section X = 
lies entirely in D_ (Fig. [2]). Let the viscous fluid (oil) in D + be sucked by a pump on the 
right infinity, with the non- viscous fluid (water) coming to D_ from the left infinity, then 
the interface T moves to the right, with normal velocity V n = V n (Z) at any point Z G T 
being given by 

V n (Z) = -d n Q(Z). (1.1) 

The potential function <3> is proportional to the pressure in the viscous fluid. It is a 
harmonic function in D_ equal to on V (zero surface tension) with the asymptotic 
behaviour $ ~ — |X as X — > +oo. The velocity field everywhere in the viscous fluid is 
given by the Darcy law V = — V$. 

The simply-connected case (a single interface) allows for an effective application of 
the conformal mapping technique (see, e.g., [21 US])- In the complex coordinates Z, Z one 
may describe the growth process in terms of a time dependent conformal map Z(W,t) 
from a fixed domain of a simple form in the "mathematical" IV-plane, say the half-strip, 
onto the evolving oil domain in the "physical" Z-plane. The interface itself is the image 
of the segment [0, 2m] of the imaginary axis. The Hele-Shaw dynamics is then translated 
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x=o 

Figure 2: Laplacian growth in a channel with periodic boundary conditions (a cylinder). 

to a nonlinear partial differential equation for the function Z(W,t), referred to as the 
Laplacian growth equation [TT] : 



21m(d a Z(ia)d t Z(ia)) = R, ere[0,27r]. (1.2) 
The harmonic moments of the oil domain D + (Richardson's moments), 

t k = - I! e- kZ l R dXdY , k>l, 

nkR J Jd+ 

are known to be constants of motion for the LG process [18J. The complimentary set of 
moments, 



e kZ ' R dXdY, k > 1, 

TK J JD- 

are time-dependent quantities. 



v k 

7TJ 



The dToda hierarchy and the dispersionless tau-function. We follow the ap- 
proach suggested in [TT] and map our problem on the cylinder to a contour dynamics in 
the radial geometry taking place in an "auxiliary physical plane" (the z-plane) with a 
non-uniform density. Contrary to the LG in the plane, the Lax function of the dToda 
hierarchy for the interface dynamics on the cylinder is not the conformal map Z(W) itself 
but the exponential function z = e z ^ w " R . The hierarchical times are Richardson's mo- 
ments tk- The LG equation (II. 2p plays the role of the "string equation" which uniquely 
characterizes the solution to the whole hierarchy. We construct the dispersionless tau- 
function for this solution, 

Fo{R', to, {tk}k>i, \tk}k>i), (1-3) 

which is a function of the harmonic moments tk, the variable to related to the area 
of the growing domain and depends on R as a parameter. The function Fo obeys the 
dispersionless Hirota equations for the dToda hierarchy 



Z2 -j e D(z 1 )D(z 2 )F _ z ^ e -dt D( Z1 )F _ Z2e -dt D(z 2 )F ^ /-^ 
Z 2 (l- e -D(zi)D(z 2 )F \ _ e dt (dt +D( Z1 )+D(z2))F ^ /-^gN 



1*1 

Z\'t 

z ^ z 

where D(z) = — — dt k , D(z) = ~r~^tk- ^ contains all information about the LG 

fc>i ^ k>i k 

process and the conformal maps in the sense that it allows one to find the complimentary 
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moments and the (inverse) conformal map W(Z) by the formulas 

dF Z 1?F, <?F 

w(z ^r-2 ig-g — SA' (L6) 

The dToda hierarchy is an example of the universal Whitham hierarchy introduced 
in [12]. It is a multi-dimensional extension of the hierarchies of hydro dynamic type [30J. 
The solutions can be parametrized [13] by canonical transformations in a two-dimensional 
phase space in such a way that any solution corresponds to a canonical pair of functions 
(the "twistor data"). In fact, this is equivalent to the characterization of the solutions 
via string equations. 

The connection with Hurwitz numbers. We also show that the function Fq is 
closely connected with the genus zero part of the generating function for the double 
Hurwitz numbers that count connected coverings of the sphere. The precise conneciton 
is as follows: 

n = g| + g ^"J^ <A»)+m-2)< H "^>^ « n n (!' 7 > 

Here /J, /2 are partitions of d = |/i| = \p,\ into £(fi) parts /ii > ji2 > • • • > fJ>£(p) > 
(respectively, into £(jl) parts /ii > J12 > • • • > > 0), H^fM, ft) * s the properly 

weighted number of topologically non-equivalent coverings / : CP 1 — > CP 1 of degree 
d having ramification points at and 00 of the types ft and ft respectively and I = 
£(ft) + £(ft) — 2 simple ramification points. The numbers i?^(/i, ft) are called the double 
Hurwitz numbers [2T]. In the series (11.71) only Hurwitz numbers corresponding to the 
genus-zero coverings enter. We see that they are basically the Taylor series coefficients 
of the dispersionless tau-function (11.31) . 



2 Lax equations for the dToda hierarchy and associ- 
ated growth processes 

2.1 Dispersionless Lax equations 

Following the approach of [llj, we start with the Lax equations for the dToda hierar- 
chy with certain reality conditions imposed. The main object is the Lax function z(w) 
represented as a Laurent series of the form 

z ( w ) =rw + a + — + — + ... (2.1) 

w w z 

The leading coefficient r is assumed to be real while all other coefficients tij are in general 
complex numbers. All the coefficients depend on deformation parameters (or "times" ) to 
(a real number) and t\, t%, £3, . . . (complex numbers) in accordance with the Lax equations 
in the Sato form 

^ = {A k (w), z(w)} , ^ = -{A^w' 1 ), z(w)}, (2.2) 
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where for any two functions of w, to the Poisson bracket is 



r , A df dg df dg 
a log w dt oto a log w 



Here and below the bar denotes the complex conjugation and f(w) = f(w). The reality 
condition implies that the second half of the Lax equations (with t^-derivatives) is ob- 
tained from the first one by complex conjugation provided w belongs to the unit circle. 
The generators of the flows are 

A k (w) = (z k (w)) + , A (w) = hgw. (2.4) 
For the dToda hierarchy, the (. . .) + -operation is 

(z») + := + 5 (*»)„■ 

Hereafter, (. . means taking the terms of the Laurent series with degrees belonging to 
the subset S G Z (in particular, (. . .) is the free term). Note that at k = equations 
(12. 2p become tautological identities. The second Lax function of the dToda hierarchy 
is z(w^ x ). The reality conditions (i.e. the requirement that its coefficients are complex 
conjugate to those of the z(w)) imply that it obeys the same Lax equations. 

By purely algebraic manipulations, one can show [TJ1 H3] that the Lax equations (12.21) 
with A k given by (12. 4p are equivalent to 

d t .A k (w) - d^Ajiw) + {Akiw),^^)} = 0, 

(2.5) 

d t .A k {w- 1 ) + dt.Mw) + {A^w-^A^w)} = 0, 



which is a dispersionless version of the "zero curvature" representation. In their turn, 
equations (12. 5p imply that the Lax equations (12.21) are compatible with each other, i.e., 
dt n {d tm z{w)) = d tm (d tn z(w)) and dt n {d tm z{w)) = d tm {d- tn z{w)) for all n,m. This means 
that the vector fields dt k , di n commute and these symbols can be understood as partial 
derivatives. 

Let w(z) be the inverse function to the Lax function z(w). In terms of the inverse 
function, the evolution equations (12. 2p acquire a simpler form: 

d\og w(z) _ dAk d\ogw(z) = _dAk_ (2 6) 

dt k dt ' dt k dt 

Here Ak = A k (w(z)), A k = A k (l/w(z)) are regarded as functions of z, and the derivatives 
are taken at fixed z. Treating A k s as functions of z, one can rewrite equations (12.51) in 
the form similar to (12. 6p : 

d_A L = d_Ak dAj_ _ _dAk (2 ?) 

dt k dtj ' dtk dtj 

Note that at j = this system coincides with (12. 6p . 
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It follows from the construction of the A^s that the expansion of A] e (w(z)) in a 
Laurent series in z is of the form Ak = z k + 0(1). More precisely, these Laurent series 
are 

A (w(z)) = logw(z) = logz - ~ d to v - Y, z~\ (2-8) 

fc>i 

A>(*)) = *>--±fy Vo -£%^-\ J>1. ( 2 -9) 
z fc>i K 

^•(^- 1 W) = ^^o + E^^ ! J>h (2.10) 
/ fc>i K 

where Vk are functions of the times such that d tj Vk = d tk Vj, d t;j Vk = dt k Vj. 

We are especially interested in the class of solutions such that z(w), for all tk in an 
open set of the space of parameters, is a univalent function in a neighborhood of infinity 
including the exterior of the unit circle. This means that z(w\) = z(w2) if and only if 
Wi = W2- From now on, we assume that z(w) belongs to this class. In this case z(w) is 
a conformal map from the exterior of the unit circle to a domain in the complex plane 
containing infinity while z(u> _1 ) is a conformal map from the interior of the unit circle 
to the complex conjugate domain. For technical reasons it is convenient to assume that 
the origin of the z-plane lies outside this domain. 



2.2 General solution to the Lax equations 

A general solution to the system of differential equations (12. 2p is available in an implicit 
form [13]. To present it, we need an extended version of the Lax formalism. 

By the definition of the Poisson bracket, logw and to form a canonical pair: 

{logw, t } = 1. 

The evolution according to the Lax equations can be regarded as a tfc-dependent canonical 
transformation from the pair (log w, t ) to another canonical pair whose first member is 
log z(w). It is quite natural to introduce the second member which we denote by M. 
Depending on the situation, we treat it either as a function of z and to or as a function of w 
and to through the composition M = M(z(w, t ),to) (it also depends on the deformation 
parameters tk). To find what is M, we note that the condition {log z, M} = 1 can 
be identically rewritten as d to M(z) = zd z logw(z,t ). It determines M up to a term 
depending only on z. The latter can be fixed by the requirement that M obeys the same 
Lax equations (12. 2\i . To wit, equation d tk M = {Ak, M} (where the derivatives are taken 
at constant w) is equivalent to 

d tk M(z) = wd w A k d t0 M(z) = zd z A k . 

Taking into account (12.81) . (12.91) . we can write 

M(z) = Y,kt k z k + to + Y,v k z- k . (2.11) 
fc>i fc>i 

This formal Laurent series represents an analytic function if the domains of analyticity 
for the functions represented by the series kt k z k and VkZ~ k overlap. The function 

k>l k>l 
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M is the dispersionless ( "quasiclassical" ) limit of the Orlov-Shulman operator [29]. Its 
geometric meaning depends on a particular solution. In a similar way, one can construct 
the conjugate Orlov-Shulman function, M(z), such that the transformation (log w, to) — > 
(log z~ 1 (w~ 1 ), M(z(w~~ 1 )) is canonical and M obeys the same Lax equations. The Lax 
equations imply that the composition of the canonical transformations 

(log z,M) ->■ (log in, t ) -»> (log z~\M) 

does not depend on i.e., it is an integral of motion. Moreover, any t^-independent 
canonical transformation (log z,M) — > (log2 _1 ,M) between the Laurent series of the 
form prescribed above generates a solution to the dToda hierarchy. See [13] for the 
detailed proof. 

More precisely, let (log f(w, to), g(w, to)) be a canonical pair: {log/, g} = 1. Suppose 
that the functions z, z, M, M are represented by Laurent series of the form given above 
and are connected by the functional relations 

l/ziw- 1 ) = f (z(w), M{z(w))) , M(z(w- 1 )) = g {z(w), M{z{w))) . (2.12) 

Then the function z(w) obeys all the Lax equations and its coefficients (as functions of 
t/s's) obey the dToda hierarchy. Conversely, any solution of the dToda hierarchy admits 
a representation of this form with some (/, g)-pair |13j . 

This construction can be made more explicit by introducing the generating function 
of the canonical transformation (logw;,to) — > (log/, g). An important class of solutions 
corresponds to the canonical transformations (logz, M) — > (logz -1 , M) defined by means 
of a generating function U(z, z) pj as follows: 

M = zd z U(z } z), M = zd- z U(z, z). (2.13) 

Here U(z, z) is an arbitrary different iable real-valued function of z, z. This form of the 
canonical transformation implies that the functions z(w) and z(w~ l ) are algebraically 
independent. These are solutions of generic type and we call them non-degenerate. They 
describe conformal maps of 2D domains with smooth boundaries. For non-degenerate 
solutions the "string equation" 

{z(w), ziw- 1 )} = 1 (2.14) 

U z g{z{w),z(w l )) 

where U Z z(z, z) = d z d z U(z, z) holds true. It is obtained by plugging M from ( 12 . 13[) into 
the canonical relation {z, M} = z. 



2.3 Contour dynamics 

Lax equations ( 12. 2 p can be treated as equations of a contour dynamics [UJ . The contour 
is the image of the unit circle, i.e., z(e ie ), < 9 < 2tt (Fig. [3]). We call it the Lax 
contour and denote it by 7. It depends on the deformation parameters. We assume 
that 7 is a non-self- intersecting curve encircling the origin. The function z(w) provides 
a time-dependent conformal map from the exterior of the unit circle onto the exterior of 
the Lax contour. 
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\w\=l 

Figure 3: The Lax contour. 

To derive equations of motion for the Lax contour, we need a general kinematic 
relation. Let (x(a, t), y(a, t)) be any parameterizations of a moving closed contour in the 
plane, then the normal velocity of the contour points is V n = ^ (d a yd t x — d a xd t y) , or, 
in the complex notation, 

V n = — (d a zd t z - d a zd t z) , (2.15) 



where dl = ^J(dx) 2 + (dy) 2 is the line element. The normal velocity V n is positive if it is 
directed to the exterior of the contour. 

Applying this formula to the Lax contour z(e ) with the specific parameterization 
a = 9 and t = to with all other t^s fixed, we get the normal velocity of the Lax contour 
7 at the points z(w), \w\ — 1: 

{z(w), zjw- 1 )} 

Vn= 2^)1 ' ( } 

Here z'(w) = d w z(w) and the Poisson bracket in the numerator is given by (12.31) . Eq. 
( I2.16P together with the string equation (I2.14p states that the normal velocity of the Lax 
contour at the point z e 7 is equal to 

XT ( X 1^(^)1 . v 

KW = 2^f7(^i)' ZGT (2 - 17) 

Eqs. (12.111) . (I2.13P allow us to express the deformation parameters through the ge- 
ometry of the moving contour: 

t k = -^—£ z- k - l (w)M{z(w))dz(w) = / z- k d z U dz , k > 1, (2.18) 
2mk J\w\=i 2mk j 7 

t = — <f M(z(w))dlogz{w) = — J> d z Udz. (2.19) 
2tu 2wi J-y 

We stress that £i, £2, • • ■ are kept constant, so they are integrals of motion for the contour 
dynamics (J27T7J). Let D be the compact domain bounded by the Lax contour and D c = 
C \ D its complement, then Green's formula implies that 

t k = --—[[ z~ k AUd 2 z, k>l, 



Auk JJd' 

(2.20) 

to = -7- II AUd 2 z , 



47r JJd 



where A = 4d z d z is the Laplace operator and d 2 z = dxdy. According to our assumption, 
the domain D contains the origin, so the integrals (I2.20p with positive k are well-defined 
(for small k a regularization at infinity is required). These formulas show that the t^'s 
are harmonic moments with the density function AU. The double integral representation 
of to implies that the density function is integrable everywhere in D. If AU is singular 
at some point (say, at z = 0), then one still may give sense to this double integral by 
introducing a cut-off (see the next section). The coefficients Vk in ( 12. lip have similar 
integral representations: 

v k = — I z k - 1 {w)M{z{w))dz{w) = — l z k d z Udz = — ff z k AUd 2 z. (2.21) 
2ni J\w\=i 2ni J-y AttJJd 

One can also define the logarithmic moment t> : 

v = ]- ff \og\z\ 2 AU{z,z)d 2 z. (2.22) 



47T JJd 

Again, this integral representation implies that AU is integrable everywhere in D, oth- 
erwise a cut-off is required. Similar to to, the moment vq is real. In the important case 
when U(z,z) depends only on zz, i.e., zd z U = zd z U, the logarithmic moment vq can be 
equivalently represented as the contour integral 



■J- (( (log Izfd.U-z^u) dz. (2.23) 



Let S(z) be the analytic continuation of the function d z U(z, z) away from the contour 
7, then S(z) = S+(z)+S-.(z), where S±(z) are analytic functions in D and D c respectively. 
They are given by the following integrals of the Cauchy type: 



(2.24) 



i rffiM^^ i, 

2m J-y C-z ^ 

Zm J-y z — c, z k ^ x 

From these formulas it follows that M(z) = zS(z). 



2.4 The dispersionless tau-function 

As is shown in [9], there exists a real-valued function F = F (t , {t^}, {ik}) such that 

dF = v dt + ^{v k dt k + v k dt k ), 

k>l 

i.e., 

dF 

v k = ^, k>0. (2.25) 
at k 

It is called the dispersionless tau-function and admits the following representation as a 
double integral over the domain D: 

F ° = ~ it 2 ' //d IId Af/(Z ' l) bg \ Z ~ X ~ C 1 Af/(C ' ^ d2zd2C (2 ' 26) 



10 



The dispersionless tau-function can be also written as 

2F = - — ff UAUd 2 z + t v + ^2(t k v k + t k Vk) ■ (2.27) 
4tt JJd 

It satisfies the dispersionless Hirota equations (ll.4j) . (II. 5ft . The conformal map w(z), 
inverse to the z(w), can be expressed through the dispersionless tau-function as follows: 

— k 

w ( z ) = ze -i^F -d t0 D(z)F ^ D ^ : =J2^- d tk . (2.28) 

k>i k 

(see [9] for details). 

Remark. As is argued in [II], F is the free energy of the model of N x N normal 
random matrices with the potential 2TZe J2k ^kZ k — U(z, z) in the iV — > oo limit. In terms 
of the eigenvalues, the partition function is given by the following iV-fold integral over 
the complex plane: 

r N ({ tj }, ft}) = ± . J] \*m-zn\ 2 II e-*^ + »E^M^) Sz r (2.29) 

iV . m<n j =1 

It is known that is, for any U(z, z), the tau-function of the 2D Toda lattice hierarchy. 
Under certain assumptions about the potential U(z, z), in the limit h 0, N — >■ oo such 
that t — is fixed, the function 

F = F (to,{ti},{fi}) = hm(ft 2 log r^}, {£-•})) (2.30) 

is the dispersionless tau-function given by (I2.26p . A matrix model representation of the 
form ( 12.291) for the Hurwitz tau-function with the potential U(z,z) oc (log(zz)) 2 (see 
below) was recently suggested in [25] . 

2.5 Example: Laplacian growth in radial geometry 

The simplest but important example is U(z, z) = zz which corresponds to the canonical 
transformation z = z _1 M, M = M (i.e., M = M = zz). The string equation (I2.14p 
acquires the form 

{z(w), ziw- 1 )} = 1. (2.31) 
In this case the normal velocity is given by 

V n (z) = ^\w'(z)\, 2 6 7 (2.32) 

Note that |w/( z )l is equal to the normal derivative <9 n log |w(z)| of the solution to the 
Laplace equation with a source at infinity and the Dirichlet boundary condition on the 
contour. Hence (12.321) is identical to the Darcy law for the dynamics of interface between 
viscous and non- viscous fluids confined in the radial Hele-Shaw cell, assuming that there 
is no surface tension at the interface. In this way we obtain the exterior LG problem in 
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which the viscous fluid occupies the non-compact exterior domain with a source or sink 
at infinity. Formulas (l2TT8j) or flX20J) state that 

t k = TT-r <f> z~ k zdz =- — ff z~ k d 2 z 
Amk J-y irk JJd c 

are harmonic moments of the exterior of the contour 7. Their conservation in the Lapla- 
cian growth dynamics was established by S.Richardson [18]. Eq. (12. 18ft states that the 
time variable t = t should be identified with area (divided by it) of the compact interior 
domain encircled by 7. The function S(z) is the Schwarz function of the contour 7 [31j. 
The dispersionless tau-function obeys the quasi-homogeneity condition [7]: 

4F == -t\ + 2t v + E( 2 - k )(^k + i k v k ) (2.33) 

fc>i 

with Vk = d tk F . Also the relation ^ ktkVk = kikVk holds. 

k>l k>l 

A more general example is U(z,z) = (zz) 1 ^ with N E Z. As it was pointed out 
in [5], the corresponding solutions describe the LG of Z^-symmetric exterior domains 
(i.e., symmetric under rotations through the angle 2tt/N), or, what is equivalent, the 
LG in a cone. At negative integer values of N one obtains the interior LG problems. In 
particular, N = —1 corresponds to the Hele-Shaw evolution of a compact domain with a 
point-like source or sink inside it. 

In the next section we consider another important example which corresponds to 
the Laplacian growth in an infinite channel with periodic boundary conditions in the 
transverse direction (i.e., in an infinite cylinder). Formally it is as a limiting case iV — > 
00 of the ZAr-symmetric solutions discussed above but the limit is rather tricky. An 
independent approach is suggested below. 

3 Laplacian growth in channel geometry 
3.1 The moving boundary value problem 

As is described in the introduction, the LG problem in the infinite cylinder (Fig. is 
translated to the following moving boundary value problem: 

' A$(Z) = in D+ 
$(Z + 2iriR) = $(Z) 

(3.1) 

$(z) = 0, z e r 

§{Z) = - \lleZ + . . . as TleZ^+oo. 

The last condition means that far away to the right from the interface the visocous fluid 
moves with constant velocity V = (^,0) (in our units time has dimension of length, so 
the velocity is dimensionless) . 
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Let W(Z) be the conformal map from D + to the right half of the mathematical W- 
plane factorized over shifts by 2iri such that W(Z + 2niR) = W(Z) + 2iri with the 
expansion as TZeZ — > +00 of the form 

W(Z) = Z/R + Y. Cke~ kZ/R . (3.2) 

fe>0 

Then the solution to the boundary value problem ( 13. ip is given by $(Z) = -^TZeW(Z). 
Since d n 1ZeW(Z) = \W'(Z)\ for all zeT, the normal velocity is 

V n (Z) = *\W'(Z)\. (3.3) 

The inverse conformal map, Z(W), is given by the series of the form 

Z{W) = RW + J2u k e~ kW . (3.4) 

fc>0 

The Laplacian growth equation (II. 2p for Z(W) can be derived with the help of the same 
kinematic identity (I2.15P applied to the contour in the Z-plane with o = —iW, then 
da/dl = \W'(Z)\, and the identity combined with the Darcy law yields 

dZ{W) dZ(-W) _ dZ{W) dZ(-W) 

~dW dt dt dW~ ~ ( } 

or, in terms of the Poisson bracket, {Z(W), Z(—W)}w,t = R- 

A remark on the LG problem in an infinite channel with rigid walls is in order. In 
this case instead of periodicity condition $(Z + 2niR) = &(Z) one should impose the 
no-flux conditions <9y$(Z) = dy^iZ + 27riR) = on the walls Z = X (the real line) and 
Z = X + 27iiR (the line Y = 2ttR). In particular, this implies that the tangent lines to 
the interface at the endpoints are orthogonal to the walls. This problem can be formally 
reduced to the problem in a cylinder of radius 2R with the additional Z2 reflection 
symmetry Y — > — Y which in the complex coordinates is the complex conjugation Z — ^ Z . 
This symmetry implies that the coefficients in (13. 4p should be real. 



3.2 A growth problem in the auxiliary physical plane 

Our strategy will be to map this problem to the auxiliary physical plane (the z-plane), 
where it becomes a radial growth problem of the type discussed above (Fig. H]). This is 
achieved by the conformal transformation 

z = r e z / R or Z = R\og(z/r ) } 

where r is some (time-independent) scale in the z-plane. The conformal maps Z(W) 
and z(w) are thus related by the formulas 

z { w ) = roe ^o^)/R ? w ( z) = e w(m og{z /r )) ^ 
Z{W) = Rlog(^-t) , W(Z) = logw(r e z/R ). 
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Figure 4: The auxiliary and real "physical" and "mathematical" planes. 

The contour T is mapped to the Lax contour 7, and their normal velocities, V n = 
and Vjf\ are connected by the formula 



dZ 



dz 



R 



(3.6) 



The domains D + and D + in the physical plane are mapped to D c and D respectively in 
the auxiliary physical plane. The X = section of the cylinder is mapped to the circle 
\z\ = To, so our assumption means that the domain D contains the disk \z\ < tq. 

According to the general method outlined in Section 2.2, let us take the following 
generating function of the canonical transformation: 



U(z,z) 



R 



log 



zz 
rl 



We consider it as a function on the auxiliary physical plane. Then 



— zz 
M = M = zd z U = zd- z U = R\og^r 



or, equivalently, 
and 



zz = r e 



d z d- z U 



2^M/R 

R_ 

zz 



(3.7) 

(3.8) 

(3.9) 
(3.10) 
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The normal velocity of the Lax contour in the z-plane is 



Vi*\z) = ^\v/(z)\, z€ 7 - (3.11) 
R 

Using (13. 6p and = 7—7 \W'(Z)\, we can find the normal velocity V^ z ' in the corre- 

z 

sponding point of the contour T: 

V^ Z \Z) = ^V^(z) = *\W'(Z)\ (3.12) 

which coincides with ( 13. 31) . The string equation ( 12.141) in the auxiliary physical plane 
reads 

{z(w), ziw- 1 )} = R~ 1 z{w)z{w~ 1 ). (3.13) 

After the change of variables z{w) = roe z( * w ^ R , z(w _1 ) = roe z( -~ w ^ R it becomes the 
Laplacian growth equation (13. 5p in the physical plane. This proves the isomorphism be- 
tween the radial growth problem in the z-plane and the physical problem on the cylinder. 

It remains to identify the time variables (deformation parameters) with moments 
of D + . The time t is identified with the t - var i a ble. In the case at hand the density 
AU(z, z) = R/\z\ 2 is non-integrable at z = 0, so some modifications in formulas from the 
previous section are necessary. Let B(r ) be the disk of radius r centered at the origin, 
then in the z-plane we have 
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R r , / \z k dz R ff d z 
t = — f]og{^)— = - — . 3.14 

Til J j v r ' Z 7T J JD\B(r ) ZZ 

Equivalently, in the Z-plane these integrals are represented as 



where Area(D^) is the area of the domain bounded by the curve T and the section 
X = 0. Note that D^ 0) is the image of B(ro) under the map from the auxiliary physical 
plane. The higher times (integrals of motion for the Laplacian growth) are 



tk 



JL. I z -* Xo JA\ d A = J±L[ e- kZ ' R XdZ = - 44 II e~ kZ / R d 2 Z. (3.16) 
mk J-y k tq j z ixikR Jr nkR JJd + 

The complimentary moments (dynamical variables for the Laplacian growth) are 

Vk = 1 I ^logf^)— = lL f e kZ/R XdZ = J± ff e kZ/R d 2 Z (317) 

iri J j \ r ' z iriR Jr ttR JJd- 

For k > 1 no regularization is required and, as it is clearly seen from the contour integral 
formulas in the z-plane, the t^s and Vk's do not depend on r . However, the integral 
(I2.22p for the logarithmic moment Vq diverges. One should cut off the integral at \z\ — ro' 

vo = - II d 2 z = 2t log r + J- I Xd 2 Z . (3.18) 

7T JJo\B(r ) ZZ TcR 2 JD { °> 
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An important relation between the moments can be derived by calculation of the 
integral I = - — : f U{z,z)z~ l dz with U given by (13. 7p in two different ways. First, by 



the Stokes formula, 



zdM — = — 



d 2 z 

log \z/r \ 2 — = v - 2t logr . 



7TJJD\B(r ) ZZ 7T J JD\B(r ) |2| 

On the other hand, comparing (13.71) and (13.81) . we see that for z G 7 it holds 

M 2 (z) 



U(z,z) 



2R 



(3.19) 



so 



— Im 2 ^ = 



1 /" 

-TBfiYl kt k z k +t + Y. v k z 

%tiJ ~1 k>l k>l 



-k\ 2dz *0 



1 



2# 



5^ H fc i; fc . 



fc>l 



Equating the results, we obtain the identity 



Rv o = -f + 2Rt o log r + J2 ktkVk- 
1 fe>i 



(3.20) 



3.3 The tau-function 



The dispersionless tau-function is given by the double integral in the z-plane over D ( I2.26L 
However, in our case the integral diverges at small \C\ and one should introduce a 
cut-off: 

F^-%\\ tf loglz-'-C-'l^ 

J J D\B(ro) J J D\B(ro) 



7T^ 



The same cut-off should be introduced in (12. 27ft : 



7T JJD\B(r ) 

In the physical plane these integrals are written as 



^0 



7T 





e -Z/R_ e -Z'/R 







2Fn 



irR 2 77d (0) 



d 2 Zd 2 Z'-t 2 logr 



X 2 d 2 Z + t Q v + ^2(t k v k + t k v k ). 

k>\ 



(3.21) 



2R 2 r r d 2 z 
2F = // [\ og (\ z \/r )} 2 — + t v + J2( t kVk + i k Vk). (3.22 

71 JJD\B(r ) \Z\ 2 ^ 



(3.23) 
(3.24) 



The general formulas v k = dF /dt k f!2.25j) hold in this case as well, so equation (I3.20j) is 
equivalent to the following relation for the first order derivatives of the function F : 



dF tl , 1 ^ , dF 

_ = ^ +2tologro + _p tt _ 



(3.25) 



The integral in 



can be simplified using the identity 



— // UAUd 2 z = — I Ud z Udz-- [[ \d z U\ 2 d 2 z 
4tt JJd 2m J do n J Jo 
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valid for any domain D in the z-plane. We notice that in our case l^f/j 2 = 2Ud z d 2 U, 
hence 

— ff UAUd 2 z = — l Ud z Udz = — l U{z, z)S(z)dz. 
An JJd 6ni j 7 6ni 7 7 

Using (I3.19p . we can rewrite the r.h.s. entirely in terms of M(z), so 

— // XJ MJd 2 z = ^— I M\z)— . 
An JJd \2niR 7 7 v z 

Plugging here the series (12.111) for M(z) and extracting the residues, we find: 

-1 If UAUd 2 z = J| + ^ £ kt k v k + ^ £ ( klt ktiv k+l + (k + l)t m v kVl ) . (3.26) 
An JJd QR R fe > 1 2R k ^ > 

(Note that the expression in the right hand side must be real although it can not be 
directly seen from its form.) Therefore, equation (I3.22p takes the form 

t 3 t 

k>l bK H k>l 

(3.27) 

~^~5 [klt k tiv k+ i + (k + l)t m v k vi). 
2K k,i>i 

It can be further simplified using relation ( 13 .20j) : 

t 3 1 1 

F o = ^ + to lo S r o + ~ ^2(hVk + hv k ) - — Yl, (kltktiVk+i + (k + l)t m v k vi\ (3.28) 
bK 1 k>\ k,i>i 

Let us examine how this function depends on r and R. As it follows from (I3.16p . 
f)3.17p . the moments t k , v k with k > 1 do not depend on r . The dependence on r comes 
from the cut-off at small distances of the formally divergent integrals for tp, Vq and the 
integral in ( 13.22[) . It is not difficult to see that 



d (l 



(log | 


A 






z 


2 



d 2 zj = -2 fc+1 (log r f. 

In particular, dt /dlogr = —2R, dv /dlogr = —AR\ogr . The full derivative of (I3.22p 
is then easily calculated to be 

dF 

-ARt \ogr . 



d\ogr 

4- t 2 lno - r n I = 

d\ogr \6R 



d ( t 3 \ 
Since — — ^- + tl\ogr = —ARt \ogr , we see that the function 



F = F -^-t 2 \ogr (3.29) 



does not depend on r : 



dF dF dF dtp 8F 8F 

H — n = tt; ^t^~R — = U. 



d\ogr <91ogr dt <91ogr <91ogr dt 
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The last equality means that the function F depends on tq and t only in the combination 
t + 2R\ogr (or rge* ^). 

The derivative dF^/dR at fixed t k can be found by the general variational method 
developed in [TJ E] but in our case a simpler argument works. First we pass to the 
dimensionless times t k = t k /R, then the new times t k are .R- independent. As is seen from 
equation (|3.2ip . F is of the form F = R 2 F , where F is R- independent. Therefore, we 
can write 

F (R, {t k }) = F (R, {Rh}) = R 2 F = R 2 F (1, {t k }). 

Next, taking the total /^-derivative of the identity R~ 2 F (R, {Rt k }) = F (l, {tfc}), we 
find: 

-2iT 3 F + R- 2 d R F + R- 2 (t d t0 F + 2fte]T i k d tk F ) = 

k>l 

or 

2F = Rd R F + t d t0 F + (tkdt k F + i k d- tk F ) , (3.30) 

fc>i 

where the partial derivative Or is taken at fixed t k and the derivatives dt k are taken at 
fixed R. Comparing with (13.271) and taking into account that v k = d tk F , we conclude 
that 

OF t^ 1 
- R 2 d R F = = JL + to J2 kt k v k + - ]T (klt k t lVk+l + {k + l)t M v k vt) . (3.31) 

OK " fc>i 1 k,i>i 

Using the notation /3 = 1/R, as in [27], we rewrite this equality in the form 

dpF = | + toE kt k d tk F + \ E (klt k tid th+l F + (k + l)t m d tk F d tl F ) (3.32) 

" k>l ^ k,l>l 

which agrees with the formula given in |27J modulo simple redefinition of the tau-function 
(more details are in the next section). The term with the double sum comes from the 
action of the cut-and-join operator in the limit of zero dispersion (the genus-zero term 
in the topological expansion). 

Passing in general formula f!2.28j) for the conformal map w(z) to the variables W, Z, 
we obtain 

1 —k 

W(Z) = Z/R + logro - - d 2 Q F - £ ^- e~* z l R d tk d k) F . (3.33) 

A k>i K 

Here we recoginize equation (II .6p from the introduction (written there at r = 1). 

Remark. For solutions with the reflection Z2-symmetry Z — )■ Z (which describe 
Hele-Shaw flows in a channel with rigid walls) all the moments t k are constrained to be 
real numbers, i.e., t k = i k . However, the vector fields d tk and di k of the dToda flows are 
transversal to the real submanifold defined by the conditions t k = i k . This means that in 
general the dToda hierarchy can not be restricted to the class of solutions with reflection 
symmetry. 
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3.4 Example: the case of non-zero to, t\ (trochoid) 



The simplest possible case is when all moments except to are zero: to — t, t^ — at 
k > 1. It corresponds to the uniform motion of the circular section X = const of the 

t 2 

cylinder with velocity \: X(t,a) = t/2, Y(t,a) = Ra. In this case v = —7: + 2t logr 

z 6R 

and the tau-function is 

t 3 

Fo= 6^ + t ° l0gr °- 

Next in order of complexity is the case to — t, t\ ^ 0, i*. — at k > 2. As is readily 
seen from (I3.16p . if only the first N moments are non-zero, then the series H3.4[) truncates 
at the iV-th term. In particular, in our case we have: 

Z{W) = RW + u + u ie ~ w , Z(-W) = -RW + u + tne w 

in the physical plane or 

z(w) = rw e^" 1 , r = r e Uo/R , u t = u x /R 

in the auxiliary physical plane (note that the latter equation has the form of the Lambert 
curve x = ye y for z~ l and w^ 1 ). Plugging this into the LG equation (13.51) . we get two 
conditions, one real and one complex: 

' 2R-Reu - <9 4 |ui| 2 = R 

< 

Riii = U1U0 ■ 

They can be easily integrated: 

R R 

Here the real parameters < k < 1, Y come from a complex integration constant. Set 
RX 2 = u + u — t, and let A be the positive root, then the time-dependent contour is 
given by 

Z(a, t) = 1R0 + 1Y0 + ^ K ) + RX{t) e~ i(T , 
where the real A = X(t) > obeys the equation 

A 2 = K V 2+ 1 (3.34) 

In terms of the Lambert function W(x) defined by the equation x = W(x)e w ^'^ we can 
write 

A 2 = _ W (- K V/*). 

In the coordinates X, Y the contour is 



X t A 2 
R = 2R + T + Xc ° Sa > 



(3.35) 

Y-Y . 
— - — = a — Asm a. 
R 
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Figure 5: The curtate cycloid X = | cos a, Y = a — \ sin a. 




Figure 6: The cycloid X = cos a, Y = a — sin a. 



At < A < 1 it is a trochoid (curtate cycloid, Fig. EJ. The initial value Ao = A(0) 
is found as a root of the equation Aq = n 2 e x o such that Ao — > as k — > 0. From 
equation (I3.34p it follows that A > 0. At A = 1 which corresponds to the critical value 
of time t c = 2R\ogn~ 1 — R, the curtate cycloid becomes the ordinary cycloid (Fig. [6]) 
with a singularity (a cusp) at a = 0. This is an example of the well known finite time 
singularities typical for the LG with zero surface tension. 

The calculation of the integral ( 13 . 16f) yields 

hence k = ro\ti\/R. The higher moments vanish. The complimentary moments Vk are 
obtained from ( 13 .20 j) (or (13 . 18[) ) and (I3.17p . All of them are in general non-zero. We 
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need them at k — 0, 1, 2: 



1 ,„ = 2i logr„ + || + |A 2 (2-A 2 ) 

»= C a2(2 - a2) ' 



p3 

-4/o o\2\ 



The dispersionless tau-function can be found with the help of (13.271) or (13.28 p . In our 
case all sums there are finite: 

t 3 t 1 
2F = t v + hvi + t x vi - ^ - t\v x - ^ t\v 2 . 

The last term is the contribution of the cut-and-join part of ( 13.271) . Plugging the values 
of the moments V/~ calculated above, we get: 

F = H + ^logr + ^A 2 (2A 4 -9A 2 + 12), (3.36) 

where A is defined as an implicit function of to,ti,ii by equation (13.341) which in terms 
of the moments takes the form 

R 2 A 2 = rlhhe^ (3.37) 

or A 2 = -W(-| ti\ 2 r 2 .e~R / R 2 ^J. This equation allows one to find partial derivatives of the 
function A: 

O to X = ^ztz — , O t ,X = zr-^z — , Or A 



2i?(l-A 2 )' 1 2tx(l-A 2 )' n 2i? 2 (l-A 2 ) 

and directly check the formulas Vq = dF /dto, V\ = dFo/dti, as well as the first equation 
of the dToda hierarchy 

dl h F = e 9 ^. (3.38) 

Note that the i^'s, as well as their derivatives with respect to to, £i, are non-singular 
at the critical point A = A c = 1 which corresponds to the critical value of to 

4 c) =2 J Rlog(-^) - it!. 

However, their second derivatives, for example 

d 2 v 1 



dtl R(l - A 2 ) 

are singular at the critical point. 

In the limit R — > oo (f3 = 0) equation (I3.37P yields A — » with RA — > r |ti|, so from 
(I3.36P we have 

F 



q = t logr + r ti*i at t k = 0, k > 2. 
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In the general case one can show that 



F = *g log r„ + £ krfhh- (3.39) 

k>l 



Set x 2 = /3 2 r 2 t 1 t 1 e /3 '°, then the Taylor expansion of the function A 2 as x — >■ reads 

oo ; 

_ \A/( r-\ = V - 

fc! 2 3 



oo h.k-1 q 

A 2 = _ W (-x 2 ) = Y / ^ Tr x 2k = x 2 + x 4 + ^x G + ^x 8 + ... 

k=l 



Plugging this into (I3.36p . we get for the restriction of the function F (/3, r^,t , {i^}, {tk}) 
to the submanifold tj — 0, j > 2, 

F {/3, r 2 , t , h, h) =§i + t 2 logr + rge^titi + \ /3 2 rfc 2 #°t?t? 

+ \ [5 A rte m n\ti + § /3 6 r 8 e 4 ^l + . . . 

At t = we get the series 

F = rgtxt! + ± /3 2 r 4 t 2 tl + I /3 4 r 6 t^ + 1 /3 6 r 8 t 4 t 4 + . . . = r 2 £ ^(/3r ) 2 - 2 ^. 
4oo n>1 n. 



As is explained in the next section, it has the form 

r 2d a2d-2 



r 2d a 

E 'O >± tt (id id\.drd 

— — H d! 2d-2{i- ,1 J tit 

d>i y za ~~ z > 



where 



d\ 

is the number of degree d coverings CP 1 — > CP 1 with exactly 2d — 2 simple ramification 
points. 

4 The LG tau-function as the generating function 
for genus zero double Hurwitz numbers 

As it was already pointed out in the introduction, the dispersionless tau-function F for 
the LG in a channel is closely related to the genus-zero part of the generating function 
for the double Hurwitz numbers which count ramified coverings of CP 1 of genus with 
arbitrary ramification type at two marked points. In this section we outline the precise 
connection between the LG tau-function ( 13. 2 1 p and the generating function for the double 
Hurwitz numbers. 

The Hurwitz numbers count ramified coverings of CP 1 . Let / : E — > CP 1 be a 
degree d covering of the Riemann sphere CP 1 by a (connected) Riemann surface E. The 
degree d of the covering is defined as the number of sheets above a generic point of CP 1 . 
A partition \i = (fix, /i 2 , • • • , fie^)) of d is a set of positive integers \ii such that 

d = ^\ii\= \n\ and ^ > ji 2 > . . . > He^) > 0. 

i=l 
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Another notation for the partition is fi = (l m i2 m2 . . .) where m« is the number of parts 
equal to i, with £(fi) being the total number of parts of the partition fi. Ramification 
points of a given <i-fold covering are classified by partitions of d which determine the type 
of monodromy (permutation of sheets) in moving around the point. The partition corre- 
sponding to a generic point (without ramification) is (l d ). A ramification point is called 
simple if the corresponding partition is of the type (l d_2 2) that means a permutation of 
two sheets. 

Let /i, p, be two partitions of d. The double Hurwitz numbers Hd,i(^,p) [21] are 
defined as the properly weighted numbers of topologically non-equivalent coverings having 
ramification points at and oo of the types \x and p, respectively, I > simple ramification 
points Pi, . . . ,Pi G CP 1 and unramified over all points other than 0, Pi, . . . , Pi, oo. The 
genus g of E is determined by the ramification data l,fi,p with the use of the Riemann- 
Hurwitz formula 

2g-2 = l-e(ji)-£(p). (4.1) 

Let F^- H \f3, Q, t, t), where t = {ti, t 2 , . . .}, t = {ti, i 2 , . . .} are the times and (3, Q are 
parameters, be the generating function for the double Hurwitz numbers: 

ai «oo m 

F^ H \P,Q,t,t) = Y,^T.Q d E H d ^,p)H^Y[fi t i Pi . (4.2) 

l>0 d>l \n\=\fi\=d i=l 8=1 

Here we employ the notation of [27] (the parameters (3, Q are related to our R, r as 
(3 = 1/ R, Q = Tq). The sum over d combined with the sum over partitions such that 
I A* I = I A* I = ^ can be written as a sum over all partitions fi, p, with the convention that 
Hd t i(fi,fi) = unless = \p,\. In [21] Okounkov has proved that the (dispersionfull) 
tau-function 

r n (t,t) = en^^K^g^^exp^^^e^^Q.t.t)) 

solves the 2D Toda lattice hierarchy of Ueno and Takasaki. 

In order to extract the contribution of genus g surfaces the following trick is usu- 
ally applied. Let us rescale {tk}, {4} and (3 by introducing a new parameter % as 
tk — > tk/h, (3 Tif3 and consider the modified generating function F( H >(h; f3, Q, t, t) := 
h 2 F^ H \h(3,Q,t/h,t/K), then the series (14. 2p having regard to the Riemann-Hurwitz for- 
mula acquires the form of the topological expansion 

FW(h;P,Q,t,t) = J2n 29 F^ H \p,Q,t,t), (4.3) 

where 

^ = ?i M 5=- ^ +'W+*9-Q ! H «MV»»*^ # n^n^ (4-4) 
counts the connected coverings of genus g. In particular, 

r^d it i — \ tip) ^(A) 

*SS3g» n^)n'^) 
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is the generating function for the numbers of the ramified coverings CP 1 — > CP 1 . 

Let F ({3, r , t , t, t) be our LG tau-function given by (I3.2ip or (13.281) . regarded as a 
function of the independent times tk, including t , and the parameters (3,r . We claim 
that the precise relation is as follows: 



6 



+ ^logr + F W (/3,r V t0 ,t,t) 



„2 fit 



(4.6) 



(equation (11 -7p from the Introduction is obtained from here at r = 1). In particular, 



F (H) (M,t,t) 



t =0 



(4.7) 



An explicit example is given at the end of the previous secrtion. 

Let us sketch a proof of equation (14. 6p . Taking into account the "initial condition" 
(I3.39p . it is enough to show that equation (I3.32p is equivalent to the "cut-and-join" 
equation for F (H) [281 [27]. We have: 



(H) 



dto 



0Q 



dF^ H) 



dQ 



=r 2 e /3t 



dF [ n H) dF^ H) 



dF^ H) 



d/3 



+t Q u±0 



d/3 
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From (14. 6 p we conclude that 
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(4.8) 



Substituting here v from (I3.2UI) . we find dFq /dt — /3 kt k Vk- Plugging this into the 

k>\ 

first equation of (14. 8p . we arrive at 



dF™ 



dF tf, ^ , 

"afl" = T + to 1^ ktkVk + 
d/3 6 ^ 9/3 



Together with (I3.32p and the fact that dF /dtk = dF^ /dtk it gives the relation 



d?Fo {H) = \ E ( WA fc+; F (H) + (k + l)t w d tk F {H) d tl F^ 



k,l>l 



(4.9) 



which is exactly the dispersionless limit of the cut-and-join operator action given in 
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5 Conclusion 



We have seen that conformal maps of plane domains and connected genus-0 ramified 
coverings of the sphere are governed by the same "master function", Fq, which is a special 
solution of the dispersionless Toda lattice hierarchy. Its arguments (commuting flows of 
the Toda hierarchy) are harmonic moments of the domain in the former case and formal 
variables necessary for constructing the generating function of Hurwits numbers in the 
latter. The double Hurwitz numbers H^i(fi,fl) for the genus-zero coverings are basically 
the coefficients of the Taylor expansion of Fq around the point tk = 0. This suggests 
that there should exist a direct connection between conformal maps and enumerative 
algebraic geometry of ramified coverings. One may also hope that this connection will 
be helpful for effectivization of the Riemann mapping theorem in the spirit of [32J. 
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